Template banks to search for low-mass binary black holes in advanced gravitational-wave detectors. 
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Coalescing binary black holes (BBHs) are among the most likely sources for the Laser Interferometer 
Gravitational-wave Observatory (LIGO) and its international partners Virgo and KAGRA. Optimal searches 
for BBHs require accurate waveforms for the signal model and effectual template banks that cover the mass 
space of interest. We investigate the ability of the second-order post-Newtonian TaylorF2 hexagonal template 
placement metric to construct an effectual template bank, if the template waveforms used are effective one 
body waveforms tuned to numerical relativity (EOBNRv2). We find that by combining the existing TaylorF2 
placement metric with EOBNRv2 waveforms, we can construct an effectual search for BBHs with component 
masses in the range 3M© < mi,m2 < 25 M©. We also show that the (computationally less expensive) Tay- 
lorF2 post-Newtonian waveforms can be used in place of EOBNRv2 waveforms when M < 12 M@. Finally, 
we investigate the effect of modes other than the dominant I — m — 2 mode in BBH searches. We find that 
for systems with (7711/7712) < 1.5, there is no significant loss in the total possible signal-to-noise ratio due to 
neglecting modes greater than / = m = 2 in the template waveforms. For higher mass ratios, including higher 
order modes could increase the signal-to-noise ratio by as much as 8% in Advanced LIGO. Our results can be 
used to construct matched-filter in Advanced LIGO and Advanced Virgo. 

PACS numbers: 04.80.Nn, 04.25.Nx, 04.30.Db, 



I. INTRODUCTION 

Over the last decade, there has been tremendous progress 
towards the first direct detection of gravitational waves. Con- 
struction of the Advanced Laser Interferometer Gravitational- 
Wave Observatory (aLIGO) is underway, with completion 
scheduled for 2014 [1]. Similar upgrades to the French- 
Italian Virgo detector [2] have commenced and construction 
of the Japanese KAGRA detector has begun [3]. When these 
second-generation gravitational-wave detectors reach design 
sensitivity, they will increase the observable volume of the 
universe by a thousandfold or more [4], compared to the first- 
generation detectors. 

The inspiral and merger of binary black holes (BBHs) 
are expected to be an important source for detection by 
aLIGO [5]. The rate of BBH coalescences that will be ob- 
served by aLIGO at design sensitivity is estimated to be be- 
tween 0.2 yr -1 and 1000 yr _1 [6]. Accurate knowledge of 
the gravitational-wave signals generated by BBHs is crucial 
for detecting and extracting information about these sources. 
To provide such waveforms, the effective one body (EOB) 
model [7] has been calibrated to numerical simulations of 
black hole mergers [8-15]. A new EOB waveform family 
(called EOBNRv2) has been recently proposed that incor- 
porates information from several non-spinning BBH simula- 
tions, with black hole ring-down quasi-normal modes [16, 17] 
attached to provide a complete BBH waveform [15]. The 
EOBNRv2 waveform is believed to be sufficiently accurate to 
search for non-spinning BBH signals in the aLIGO sensitive 
band (10-1000 Hz). 

Past searches for BBHs [18-22] used matched-filtering [23, 
24] to search for coalescing compact binaries. These searches 
divided the BBH mass space into a low-mass region with 
M = mi + m2 ^ 25 Mq and a high-mass region with 
M > 25 M Q . In this paper, we focus attention on BBH sys- 
tems with component masses between 3 M < mi , 7772 < 
25M , which encompasses mass distribution of black hole 



candidates observed in low-mass X-ray binaries [25]. aLIGO 
will be able to detect coalescing BBH systems with compo- 
nent masses mi = m 2 = 25 M to a maximum distance of 
up to ~ 3.6 Gpc. Since we do not know a priori the masses of 
BBHs that gravitational- wave detectors will observe, searches 
use a bank of template waveforms which covers the range of 
BBH component masses of interest [26, 27]. This technique 
is sensitive to the accuracy of the waveform templates that are 
used as filters and the algorithm used to place the template 
waveforms [28]. An accurate template bank is required as in- 
put for matched filter searches in the Fouier domain [24], as 
well as newer search algorithms such as the singular value de- 
composition [29]. 

In this paper, we investigate three items of importance to 
advanced-detector BBH searches: First, we study the ac- 
curacy of template placement algorithms for BBH searches 
using EOBNRv2 waveforms. Optimal template placement 
requires a metric for creating a grid of waveforms in the 
desired region of parameter space [30], however no ana- 
lytic metric exists for the EOBNRv2 waveform. In the ab- 
sence of such a metric, we construct a template bank us- 
ing the second-order post-Newtonian hexagonal placement 
algorithm [31-34]. This metric is used to place template 
grid points for the aLIGO zero-detuning high power sensi- 
tivity curve [35] and we use EOBNRv2 waveforms at these 
points as search templates. We find that the existing algorithm 
works well for BBHs with component masses in the range 
3M < mi, m2 < 25 Mq. For a template bank constructed 
with a minimal match of 97% less than 1.5% of non-spinning 
BBH signals have a mismatch greater than 3%. We therefore 
conclude that the existing bank placement algorithm is suffi- 
ciently accurate for non-spinning BBH searches in this mass 
region. Second, we investigate the mass range in which the 
(computationally less expensive) third-and-a-half-order Tay- 
lorF2 post-Newtonian waveforms [26, 36-44] can be used 
without significant loss in event rate, and where full inspiral- 
merger-ringdown EOBNRv2 waveforms are required. We 
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construct a TaylorF2 template bank designed to lose no more 
than 3% of the matched filter signal-to-noise ratio and use the 
EOBNRv2 model as signal waveforms. We find that for non- 
spinning BBHs with M < 12 M©, the TaylorF2 search per- 
forms as expected, with a loss of no more than 10% in the 
event rate. For higher masses larger event rate losses are ob- 
served. A similar study was performed in Ref. [45] using an 
older version of the EOB model and our results are quantita- 
tively similar. We therefore recommend that this limit is used 
as the boundary between TaylorF2 and EOBNRv2 waveforms 
in Advanced LIGO searches. Finally, we investigate the ef- 
fect of modes other than the dominant I = m = 2 mode 
on BBH searches in aLIGO. The horizion distance of aLIGO 
(and hence the event rate) is computed considering only the 
dominant mode of the emitted gravitational waves, since cur- 
rent searches only filter for this mode [6] . However, the inclu- 
sion of sub-dominant modes in gravitational-wave template 
could increase the reach of aLIGO [46, 47]. If we assume 
that BBH signals are accurately modeled by the EOBNRv2 
waveform including the five leading modes, we find that for 
systems with (rai/7712) < 1.5, there is no significant loss in 
the total possible signal-to-noise ratio due to neglecting modes 
other than I = m = 2 in the template waveforms, if one uses 
a 97% minimal-match bank placed using the TaylorF2 met- 
ric. However, for higher mass ratios, including higher order 
modes could increase the signal-to-noise ratio by as much as 
8% in aLIGO. This increase in amplitude may be offset by the 
increase in false alarm rate from implementing searches which 
also include sub-dominant waveform modes in templates, so 
we encourage the investigation of such algorithms in real de- 
tector data. 

The remainder of this paper is organized as follows: In 
Sec. II we review the gravitational waveform models used in 
this study. In Sec. Ill we present the results of large-scale 
Monte Carlo signal injections to test the effectualness of the 
template banks under investigation. Finally in Sec. IV we re- 
view our findings and reccomendations for future work. 



II. WAVEFORMS 

The dynamics of a BBH system can be broadly divided 
into three regimes: (i) The early inspiral, when the separation 
between the black holes is large and their velocity is small, 
can be modeled using results from post-Newtonian (PN) the- 
ory [48]. The gravitational-wave phasing of non- spinning bi- 
naries is available up to 3.5PN order [38^4]. (ii) Accurately 
modeling the late-inspiral and merger requires the numeri- 
cal solution of the Einstein equations [49-55]. (iii) The fi- 
nal ring-down phase can modeled using a super-position of 
quasi-normal modes (QNMs) which describe the oscillations 
of the perturbed Kerr black-hole that is formed from the coa- 
lescence [16, 17]. 

Numerical simulation of BBH systems are computationally 
expensive, and results are only available for a relatively small 
number of binary systems (see e.g. [56]). The EOB model [7] 
provides a framework for computing the gravitational wave- 
forms emitted during the inspiral and merger of BBH systems. 



By attaching a QNM waveform and calibrating the model to 
numerical relativity (NR) simulations, the EOB framework 
provides for accurate modeling of complete BBH waveforms 
(EOBNR). The EOBNR waveforms can be computed at rel- 
atively low cost for arbitrary points in the waveform param- 
eter space [8-15]. In particular the EOB model has recently 
been tuned against high-accuracy numerical relativity simula- 
tions of non-spinning BBHs of mass-ratios q = {1, 2, 3, 4, 6}, 
where q = rai/7712 [15]; we refer to this as the EOB- 
NRv2 model, which we review the major features of below. 
Throughout, we set G = c = 1. 

The EOB approach maps the fully general-relativistic dy- 
namics of the two-body system to that of an effective mass 
moving in a deformed Schwarzschild spacetime [7]. The 
physical dynamics is contained in the deformed- spacetime 's 
metric coefficients, the EOB Hamiltonian [7], and the 
radiation-reaction force. In polar coordinates (r, <£), the EOB 
metric is written as 

ds 2 ff = -A(r)dt 2 + §^dr 2 + r 2 (dO 2 + sin 2 6d<D 2 ) . 

(1) 

The geodesic dynamics of the effective mass /i = mira^/M 
in the background of Eq. (1) is described by an effective 
Hamiltonian H eE [7, 57]. The EOBNRv2 model uses Pade- 
resummations of the third-order post-Newtonian Taylor ex- 
pansions of the metric coefficients A(r) and D(r), with addi- 
tional 4PN and 5PN coefficients that are calibrated [9-12, 15] 
to ensure that the dynamics agrees closely with NR simula- 
tions of comparable mass binaries. 

Gravitational waves carry energy and angular momentum 
away from the binary, and the resulting radiation-reaction 
force causes the orbits to shrink. This is related to the 
energy flux as 
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where, v = = (irMf) 1 / 3 and / is the instantaneous 

gravitational- wave frequency. The energy flux dE/dt is ob- 
tained by summing over the contribution from each term in 
the multipole expansion of the waveform, i.e. 
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1Z is the physical distance to the binary, and hi m are the multi- 
poles of the waveform when it is decomposed in spin weighted 
spherical harmonic basis according to 
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where Y l J^ are the spin weighted spherical harmonics, and 
h+ and h x are the two orthogonal gravitational wave polar- 
izations. These waveform multipoles depend on the coordi- 
nates and their conjugate momenta, and their Taylor expan- 
sions were re-summed as products of individually re-summed 
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where e is if (I + m) is even, and is 1 otherwise. This 
factorized-re-summation of the waveform multipoles ensures 
agreement with NR waveform multipoles [8-10]. The first 
factor /^m'^ i s me re-summation of the Newtonian order con- 
tribution and the second factor is the source term, given 
by the mass or the current moments of the binary in the EOB 
formalism [58, 59]. The tail term 7] m is the re-summation of 
the leading order logarithmic terms that enter into the trans- 
fer function of the near-zone multipolar waves to the far- 
zone [59]. The last term Ni m attempts to capture the non- 
circularity of the quasi-circular orbits. While calculating the 
energy flux in this study we follow exactly the prescription 
of Ref. [15], which calibrates the coefficients of the flux so 
that resulting EOB waveform multipoles reproduce their NR 
counterparts with high accuracy. 

We use the EOBNRv2 Hamiltonian and flux in the equa- 
tions of motion for the binary, given by 
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where, i(= t/M) is time in dimensionless units. Starting 
from the initial coordinate configuration calculated for a fidu- 
ciary starting gravitational wave frequency, we numerically 
integrate these equations to get the evolution of the coordi- 
nates and momenta (r(t), Q(i),p r (i),p<$>(i)) over the course 
of inspiral, until the light-ring is reached. In the EOB model, 
the light-ring is defined as the local maximum of the orbital 
frequency Vt. From the coordinate evolution, we also calcu- 
late hf m (t), which is the analytic expression for the wave- 
form multipole without the non-quasi-circular correction fac- 
tor (defined in Eq. (5b)). While generating hf m (t) from the 
dynamics, the values for the free parameters in the expressions 
for Si m and pi m , are taken from Eqn.[38a-19b] of Ref. [15], 
where they optimize these parameters to minimize the phase 
and amplitude discrepancy between the respective EOB wave- 
form multipoles and those extracted from NR simulations. 

The EOB ringdown waveform is modeled as a sum of TV 
quasi-normal-modes (QNMs) [9, 10, 12, 16] 
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where N = 8 for the model we consider. The matching time 
is the time at which the inspiral-plunge and the ring- 
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at which the amplitude of the inspiral-plunge part of hi m (i) 

peaks ^i.e.tp^ ak ^ [9, 15]. The complex frequencies of the 

modes (Jimn depend on the mass Mf and spin aj of the BH 
that is formed from the coalescence of the binary. We use the 
relations of Ref. [15], given by 



0.4333^-0.4392^, (11a) 
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Using the mass and spin of the final BH, the complex fre- 
quencies of the QNMs can be obtained from Ref. [16], where 
these were calculated using perturbation theory. The com- 
plex amplitudes Ai mn are determined by a hybrid-comb nu- 
merical matching procedure described in detail in Sec. II C of 
Ref. [15]. 

Finally, we combine the inspiral waveform multipole 
him (t) and the ringdown waveform h RB (t) to obtain the com- 
plete inspiral-merger-ringdown EOB waveform h lMR (t), 



hZ K (t) = hi m {t)Q{t\ 



match 



t)+h tlu (t)e(t-tz t n, (12) 



where Q(x) = 1 for x > 0, and otherwise. These multipoles 
are combined to give the two orthogonal polarizations of the 
gravitational waveform, h+ and h x , 



l m 



(13) 



down waveforms are attached and is chosen to be the time 



where i is the inclination angle that the binary's angular mo- 
mentum makes with the line of sight, and C is a fiduciary 
phase angle. To ensure the correctness of our results, we 
wrote independent code to implement the EOBNRv2 wave- 
form based solely on the content of Ref. [15]. We then vali- 
dated our code against the EOBNRv2 waveform algorithm in 
the LSC Algorithm Library (LAL) [60]. We find agreement 
between these two implementations, giving us confidence in 
both our results and the correctness of the LAL EOBNRv2 
code. 

Previous searches for stellar-mass BBHs with total mass 
M < 25M in LIGO and Virgo used the restricted Tay- 
lorF2 PN waveforms [26, 36, 37]. Since this waveform is 
analytically generated in the frequency domain, it has two 
computational advantages over the EOBNRv2 model: First, 
the TaylorF2 model does not require either the numerical so- 
lution of coupled ODEs or a Fourier transform to generate 
the frequency domain signal requred by a matched filter. We 
compared the speed of generating and Fourier transforming 
EOBNRv2 waveforms, to the speed of generating Taylor F2 
waveforms in the frequency domain, and found that the for- 
mer can be O(10 2 ) times slower than the latter. Second, the 
TaylorF2 model can be implmented trivially as a kernel on 
Graphics Processing Units, allowing search pipelines to lever- 
age significant speed increases due to the fast floating-point 
performance of GPU hardware. We found the generation of 
TaylorF2 waveforms using GPUs to be (9(1 4 ) times faster 
than generating and Fourier transforming EOBNRv2 wave- 
forms on CPUs. However, use of the TaylorF2 waveform 
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may result in a loss in search efficiency due to inaccuracies 
of the PN approximation for BBHs. To investigate the loss in 
search efficiency versus computational efficiency, we use the 
restricted TaylorF2 waveform described below. 

The Fourier transform of a gravitational waveform h(t) is 
defined by 



Hf) = / 



-27ri 



ft h(t)dt. 



(14) 



Using the stationary phase approximation [61], the Taylor F2 
waveform h(f) can be written directly in the frequency do- 



main as 
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(15) 



where we have kept only the leading-order amplitude terms; 
this is known as the restricted PN waveform. The amplitude 
A oc M b J Q /1Z, where M c is the chirp-mass of the binary, 
A4 C = (mi + rri2)rf^, rj = ra 17712/(7711 + m^) 2 is the 
symmetric mass ratio, and 1Z is the distance to the binary. The 
Fourier phase of the waveform at 3.5PN order is given by [24, 
26, 36, 62-65] 
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where v = (ttM/) 1 / 3 is the characteristic velocity of the bi- 
nary, and 7 is Euler's constant. The initial conditions are set 
by starting the waveform from a given gravitational- wave fre- 
quency / = /low and the waveform is terminated at the fre- 
quency of a test particle at the innermost stable circular orbit 
(ISCO) of a Schwarzschild black hole (r = 6M). 



where S n (f) is the one-sided power spectral density (PSD) 
of the detector noise. Throughout, we use the aLIGO zero- 
detuning high power noise curve as the PSD for bank place- 
ment and overlap calculations, and set / m i n = 15 Hz. The 
normalized overlap between the two waveforms is given by 



(fti|ft 2 ) 



(fti|ft 2 ) 



(18) 



In addition to the two mass-parameters of the binary, this nor- 
malized overlap is also sensitive to the relative phase of coa- 
lescence (j) c and to the difference in the time of coalescence 
between the two waveforms hi and h^tc These two param- 
eters (0c and t c ) can be analytically maximized over to get the 
maximized overlap O 

0(fti,ft 2 ) = max (fti|ft 2 (0 c ,*c)) , (19) 

<f> c , t c \ J 

which gives a measure of how "close" the two waveforms are 
in the waveform manifold. The mismatch M between the 
same two waveforms is written as, 



M(h 1 ,h 2 ) = l-0(h 1 ,h 2 ). 



(20) 



If h e a is the waveform emitted by a BBH system then the Fit- 
ting Factor of a bank of template waveforms (modeled using 
approximant X) for this waveform, is defined as the maximum 
value of maximized normalized overlaps between h e a and all 
members h£ of the bank of template waveforms [28]; i.e. 



G bank 



(21) 



This quantity simultaneously quantifies the loss in recovered 
signal-to-noise ratio (SNR) due to the discreteness of the 
bank, and the inaccuracy of the template model. The similarly 
defined quantity MM (minimal match) quantifies the loss in 
SNR due to only the discreteness of the bank as both the exact 
and the template waveform is modeled with the same wave- 
form model, i.e. 



MM = min max 0(h*,h£), 

a be bank 



(22) 



where a is any point in the space covered by the bank, and X 
is the waveform approximant. For a detection search that aims 
at less than 10% (15%) loss in event detection rate due to the 
discreteness of the bank and the inaccuracy of the waveform 
model, we require a bank of template waveforms that has TT 
above 0.965 (0.947) [45, 66, 67]. 



III. RESULTS 

To assess the effectualness of the template banks con- 
structed here, we compute the fitting factors [28] of the tem- 
plate bank, defined as follows. The frequency weighted over- 
lap between two waveforms h\ and h 2 is given by 



(fti|ft 2 ) = 2 
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u J mil 



Sn(f) 



df, (17) 



A. EOBNRv2 templates placed using TaylorF2 metric 

In this section we measure the effectualness of the second- 
order post-Newtonian hexagonal template bank placement 
metric described in Ref. [32] when used to place EOBNRv2 
waveform templates for aLIGO. The same template place- 
ment algorithm was used to place a grid of third- and- a-half or- 
der post-Newtonian order TaylorF2 waveforms for low-mass 
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FIG. 1. This figure shows the effectualness of a bank of EOBNRv2 
templates, placed using the 2PN accurate hexagonal template place- 
ment of Ref. [32], to search for a population of BBH signals sim- 
ulated with EOBNRv2 waveforms. The masses of the BBH popu- 
lation are chosen from a uniform distribution of component masses 
between 3 and 25 M . For each injection, we plot the component 
masses of the injection, and the fitting factor (J 7 J 7 ). 



BBH detection searches for initial LIGO and Virgo observa- 
tions [18-22]. We construct a template bank which has a de- 
sired minimal match of 0.97 for waveforms with component 
masses between 3M < mi, m 2 < 25M . This template 
bank contains 10, 180 template grid points in (mi, 7712) space 
for the aLIGO noise curve, compared to 373 grid points for 
the initial LIGO design noise curve. For the template wave- 
forms at each grid point, we use the EOBNRv2 waveforms, 
rather than TaylorF2 waveforms. Since the metric itself was 
derived using second-order TaylorF2 waveforms, we do not, 
a priori know if this metric is a good measure to use to place 
template banks for EOBNRv2 waveforms. 

To test the effectualness of this template bank, we perform 
a Monte-Carlo simulation over the 3M <mi,m2 < 25M 
BBH mass space to find regions where the bank placement 
algorithm leads to under-coverage. We sample 90,000 points 
uniformly distributed in individual component masses. For 
each of these points, we generate an EOBNRv2 waveform 
for the system with component masses given by the coordi- 
nates of the point. We record the TT of the template bank 
for each of the randomly generated BBH waveforms in the 
Monte-Carlo simulation. Since we use EOBNRv2 waveforms 
both to model the true BBH signals and as matched-filter tem- 
plates, any departure in fitting factor from unity is due to the 
placement of the template bank grid. 

For a bank of template waveforms constructed with a MM 
of 0.97, Fig. 1 and Fig. 2 show that the TT of the bank re- 
mains above 0.97 for ~ 98.5% of all simulated BBH sig- 
nals. Less than ~ 1.5% of signals have a minimal match of 
less than 0.97, with the smallest value over the 90,000 sam- 



FIG. 2. This figure shows a cumulative histogram of the fraction of 
the BBH signal space (on the y-axis), where the bank of EOBNRv2 
waveforms has TT less than the respective values on the x-axis. The 
EOBNRv2 bank has a fitting factor TT below the 0.97 for less than 
~ 1.5% of all simulated signals with component-masses rai,ra2 
between 3 M© and 25 M© . 



pled points being ~ 0.96. The diagonal features observed in 
Fig. 1 are due to the hexagonal bank placement algorithm and 
are related to the ellipses of constant chirp mass in Fig. 4 of 
Ref. [32]. From these results, we conclude that the existing 
template bank placement metric adequately covers the BBH 
mass space with EOBNRv2 waveform templates; it is not nec- 
essary to construct a metric specific to the EOBNRv2 model. 
aLIGO detection searches can employ the second-order post- 
Newtonian bank placement metric with the hexagonal place- 
ment algorithms [30-34] to place template banks for EOB- 
NRv2 waveforms without a significant drop in the recovered 
signal-to-noise ratio. 



B. Effectualness of TaylorF2 templates 

We next explore the efficiency of using the computation- 
ally cheaper TaylorF2 waveforms to search for a population 
of BBH signals with component masses between (3-25)M . 
The signals from this population are modeled with the full 
EOBNRv2 waveforms. We use the same template bank place- 
ment as above, however now we use the third-and-a-half PN 
order TaylorF2 model as the template waveforms. This model 
does not capture the merger and ringdown of BBH signals, as 
it is terminated at the Schwarzchild test-particle ISCO. Fur- 
thermore, it diverges from the true BBH signal in the late in- 
spiral. It is important to determine when these effects become 
important. 

We sample the (3-25) M BBH component mass space at 
25,000 points by generating an EOBNRv2 waveform to gen- 
erate the "true" signal waveform. We generate a bank of Tay- 
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FIG. 3. The fitting factor TT of a bank of Taylor F2 waveforms, 
constructed with MM = 0.97 used to search for BBH systems mod- 
eled using EOBNRv2 signals. Each point corresponds to a simulated 
signal and the color bar shows the measured TT for that signal. 



FIG. 4. This figure shows the cumulative histogram of the the fitting 
factors of a bank of TaylorF2 waveforms against EOBNRv2 wave- 
forms. The points shown here cover only the component mass region 
with total masses below 12M©. 



lorF2 template waveforms over the same region, and calculate 
its TT for each of the sample points, against the correspond- 
ing EOBNRv2 waveform. Fig. 3 shows the distribution of 
the TT obtained for the TaylorF2 bank. Clearly the TaylorF2 
bank is not effectual for the enture BBH region considered, 
with mismatches of up to 18% observed. In Fig. 4 we re- 
strict the region of the BBH space to signals with total mass 
M < 12M . We find that the TaylorF2 template bank has 
TT above ~ 0.965 for > 99.5% of the points sampled in 
this restricted region. We conclude that the TaylorF2 bank is 
effectual for BBH signals below 12M . 

The value of our limit on total-mass is in agreement with 
the previous study in Ref. [45], however this analysis used 
the EOBNRvl model [68] and an older version of the Ad- 
vanced LIGO noise curve [45]. This agreement provides con- 
fidence that this limit will be robust in aLIGO searches and 
we propose this limit as the upper cutoff for the computation- 
ally cheaper TaylorF2 search. To investigate the loss in the 
TT due to the mismatch in the template and signal waveform 
models, we also performed a Monte-Carlo simulation using a 
denser TaylorF2 bank with MM = 0.99. We found that us- 
ing this dense bank of third-and-a-half order TaylorF2 wave- 
forms, we can relax the limit on the upper mass to 16M and 
still achieve a TT above 0.965, for over 99.8% of the sig- 
nals sampled in the region. However, increasing the minimal 
match increases the size of the template bank from 10, 180 to 
27, 882 templates. This is a significant increase, compared to 
the cost of filtering with EOBNRv2 templates. 



C. Effect of sub-dominant modes 

Having established that the second-order post-Newtonian 
hexagonal template bank is effectual for placing a bank of 



EOBNRv2 templates, we now investigate the effect of ne- 
glecting sub-dominant modes in BBH searches. The sen- 
sitivity reach of the aLIGO detectors is normally computed 
assuming that the search is only sensitive to the dominant 
I = m = 2 mode of the gravitational waveform. For 
binary black hole signals, sub-dominant modes may con- 
tain significant power [47]. A search that includes these 
modes could, in principle, have an increased reach (and hence 
event rate) compared to a search that only uses the dominant 
mode. The EOBNRv2 model of Ref. [15] has been calibrated 
against higher order modes from numerical relativity simu- 
lations. We investigate the effect of ignoring these modes 
in a search by modeling the BBH signal as an EOBNRv2 
signal containing the dominant and sub-dominant multipoles: 
him = ^22, ^21, ^33, ^44, ^55 (which we call EOBNRv2HM) 
and computing the fitting factor of leading-order EOBNRv2 
templates placed using the TaylorF2 metric. 

We simulate a population of BBH signals by sampling 
100,000 systems uniformly in the mi,m2 G [3, 25] M 
component-mass space. These EOBNRv2HM signals are uni- 
formly distributed in sky-location angles and inclination and 
polarization angles, which appear in the detector's response 
function to the gravitational- wave signal [69]. The template 
bank is again placed with a desired minimal match of 0.97 and 
for each of signal waveforms, we calculate the TT against the 
entire bank of EOBNRv2 waveform templates. Fig. 5 shows 
the value of the TT of the bank of EOBNRv2 waveform tem- 
plates at the sampled mass values. As expected, the higest 
fitting factors are observed close to the equal mass line, since 
when the mass ratio is close to unity, the amplitude of the 
sub-leading waveform modes is several orders of magnitude 
smaller than the amplitude of the dominant mode. As the mass 
ratio increases, the relative amplitude of the sub-leading mul- 
tipoles increases, as illustrated by Fig. 1 of Ref. [15] and the 
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fitting factor decreases. For systems with mass-ratio q above 
1.5 (3.25), we find that the TT of the EOBNRv2 waveform 
bank is above 0.965 (0.947) over 99.9% of this restricted re- 
gion. These results demonstrate that the effect of ignoring 
sub-dominant modes does not cause a significant loss in the 
total possible signal-to-noise if the mass ratio is less than 1.5. 
Fitting factors as low as 0.92 are observed in the most asymet- 
ric region of the mass space, suggesting that a search that in- 
cludes higher order modes could achieve a non-trivial increase 
in sensitivity over leading-order mode templates. However, 
this needs be evaluated carefully: an algorithm that includes 
sub-dominant modes could have an increased false-alarm rate 
(background) over a search that includes only the leading- 
order mode. Furthermore, it has been found that binaries for 
which the lowest fitting factors are observed radiate the least 
energy and hence are those to which the search is least sen- 
sitive [47]. Therefore the overall event rate loss from a real 
BBH signal population may not be as significant as suggested 
by fitting factor calculations alone. 

IV. CONCLUSIONS 

We used the TaylorF2 second-order post-Newtonian hexag- 
onal placement algorithm of Refs. [32-34] to construct a 
template bank of EOBNRv2 waveforms with MM of 0.97. 
We calculated the fitting factor {TT) of this bank against ~ 
90, 000 simulated EOBNRv2 signals with component masses 
uniformly distributed between 3M < mi,ra2 < 25M . 
We find that the TT of the template bank is greater than 0.97 
for 98.5% of the simulated EOBNRv2 signals, assuming the 
zero-detuning high power noise spectrum for aLIGO sensitiv- 
ity [35]. We conclude that the existing placement algorithm 
is effectual for use in aLIGO BBH searches, assuming that 
EOBNRv2 is an accurate model of BBH signals in this mass 
region. We then demonstrated that use of the computationally 
cheaper third-and-a-half order TaylorF2 waveform results in 
a loss in search efficiency due to inaccuracies of the post- 
Newtonian approximation, and neglect of merger-ringdown 
for BBHs with a total mass M > 12M . However, below 
this limit the TaylorF2 model is an acceptable signal for BBH 
serches. This was done using a bank with a MM of 0.97. 
By increasing the density of the bank to 0.99MM, the limit 
on total-mass can be relaxed to 16 M , with an increase in 
computational cost due to the number of templates increas- 



ing by a factor of ~ 2.7. Finally, we investigated the loss 
in the SNR incurred by the using template banks constructed 
using only the leading order mode of EOBNRv2 waveforms. 
We found that a leading-order EOBNRv2 template bank con- 
structed with a MM of 0.97 is effectual for BBHs searches for 
which 1 < (mi/m 2 ) < 1.5 and there is no significant loss in 
potential signal-to-noise ratio for systems as high as q = 3.25. 
For such systems, using EOBNRv2HM templates is unlikely 
to give a significant increase in the range to which they can be 
detected. For BBH with (m\/m2) > 3.25, detection searches 
could possibly gain sensitivity by the use of EOBNRv2HM 
waveforms, if they can be implemented without increasing the 
false alarm rate. 

Our results suggest that a significant portion of the non- 
spinning stellar-mass BBH parameter space can be searched 
for using LIGO's existing search algorithms. For systems with 
total mass below 12M template banks of TaylorF2 can be 
used without significant loss in event rate. For higher mass 
systems, neglecting high-order modes in an EOBNRv2 search 
does not cause a substantial reduction in the maximum possi- 
ble reach of BBH searches. Finally, we note that our study 
does not consider BBH systems with BH masses higher than 
M = 25 M , or the effect of black hole component spins. 
Future work will extend this study for systems with spinning 
and/or precessing black holes and consider the effect of non- 
Gaussian transients in real detector noise. 
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